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Abstract 

We analyze the phenomenon of nonlinear stochastic resonance (SR) in noisy bistable systems 
driven by pulsed time periodic forces. The driving force contains, within each period, two pulses 
of equal constant amplitude and duration but opposite signs. Each pulse starts every half-period 
and its duration is varied. For subthreshold amplitudes, we study the dependence of the output 
signal-to-noise ratio (SNR) and the SR gain on the noise strength and the relative duration of 
the pulses. We find that the SR gains can reach values larger than unity, with maximum values 
showing a nonmonotonic dependence on the duration of the pulses. 
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In recent work , we have carried out a detailed analytical and numerical study of 
the nonlinear response of a noisy bistable system subject to a subthreshold time periodic 
driving force. The force considered in those works is such that it remains constant within 
the duration of each half period while switching its sign every half period. We have focused 
our attention on the analysis of the dependence of the output signal-to-noise ratio and the 
corresponding SR gain on the noise strength. The analytical study was based on a two-state 



approximation amenable to exact treatment We showed that, for subthreshold input 
signals of sufficiently long periods, the phenomenon of SR can be accompanied by SR gains 
larger than unity. This is a genuine characterization of nonlinearity^asJSR gains larger than 
unity are strictly forbidden within a linear response description 
results were corroborated by numerical simulations 2L 



In the last few years, Gingl and collaborators 



BOB 



as SR g 



The analytical 



have carried out analog simulations 



of noisy bistable systems subject to driving forces of the type given in Eq. (JTJ) 



F(t) 



A : < t < t c 

-A : f < t < f + t c (1) 
: otherwise 



It is convenient to introduce the parameter r = 2t c /T, measuring the fraction of a period 
during which this driving force has a nonvanishing value (the parameter r in the present 
paper corresponds exactly to what Gingl and coworkers term "duty cyle"). In j^B,]^, the 
SNR and the SR gain for subthreshold amplitude input signals with r < 0.3 were studied. 
These authors find SR gains larger than unity, and, also, that increasing the r value lowers 
the SNR gain. They rationalize their observations by noting that the input SNR increases 
as r increases, while the output SNR is less sensitive to the value of r. The case studied 
by us in corresponds to the largest possible value of the parameter r, namely, r = 1. 

It seems therefore of interest to extend our analysis to input signals with r < 1 in order to 
compare with the predictions of Gingl et al. 

Let us consider a system characterized by a single degree of freedom, x, whose dynamics 
(in dimensionless units) is governed by the Langevin equation 

x(t) = -U'[x(t) 1 t]+£(t), (2) 

where £(i) is a Gaussian white noise of zero mean with (£(t)£(s)) = 2D5(t — s), and —U'(x, t) 
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represents the force stemming from the time-dependent, archetype bistable potential 

U(x,t) = ^-^-F(t) Xl (3) 

with F(t) given by Eq. The one-time correlation function is defined as 

C(r) = ^J o T dt(x(t + r)x(t)} 00 . (4) 

It can be written exactly as the sum of two contributions: a coherent part, C co h{T), which is 
periodic in r with period T, and an incoherent part, Cj nco / l (r), which decays to for large 
values of r and reflects the correlation of the output fluctuations about its average (the noisy 
part of the output). The coherent part C coh (r) is given by 

C coh {r) = l- /" T dt(^ + r)) oo (a:(t)) 0O , (5) 



T 

and C inco h{r) is obtained from the difference of Eq. (jlj) and Eq. (jSJ). In the expressions 
above, the susbcript indicates that the averages are to be evaluated in the limit t — > oo. 

According to McNamara and Wiesenfeld ll|, the output SNR, R out , is defined in terms 
of the Fourier transform of the coherent and incoherent parts of C(r) as 

_ lim £ ^ 0+ /g+ £ £ du C(u) 

^incoh \* -) 

where H(oj) denotes the Fourier cosine transform of H(t), i.e., H(u) = 
2/tt J °° dr H{t) cos{ujt). Note that this definition of the output SNR differs by a factor 
2, stemming from the same contribution at uj = —Q, from the definitions used in earlier 
wo* B fi. The peno dlcity of the eohe reI1 t part gi ve 8 ri se to delta pea ks in the speotnun. 
Thus, the only contribution to the numerator in Eq. (jUJ) stems from the coherent part of the 
correlation function. The output SNR can then be expressed as 

Rout = %r-, (7) 
(41 

where 

Q u = l f T drC coh (r) cos(fir), (8) 
1 Jo 

and 

2 f 00 

Qi = - dTC incoh (r) cos(fir). (9) 

7T JO 
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A nonmonotonic behavior of the SNR with the strength of the noise is a signature of the 
phenomenon of SR. 

The signal-to-noise ratio for an input signal F(t) + is given by 

R _ \Ul + 9l) nn x 
JMnp — 2jj ■ \ LK} ) 

7T 

where 

2A 

/i = — sinfit c (11) 

7T 

and 

2A 

g x = — (l-cosnt c ) (12) 

71 

t is then clear that for fixed values of A and D, R inp increases with r as pointed out in Ref. 
8]. Also, for given A and r, R inp decreases as D increases. 

The SR gain, G, is defined as the ratio of the SNR of the output over that of the input; 
namely, 

G=|K (13) 



stochastic resonators with static nonlinearities driven by square pulses 



SR gain values larger than I have been obtained in driven nondvnaniical systems in 



121 ]. or in noisy 



bistable systems driven by superthreshold input sinusoidal frequencies 13]. The existence 
of SR gains with values larger than 1 indicates a truly nonlinear SR. 

Although the two-state approximation introduced in Ref. p| can, in principle, be ex- 
tended to analyze systems driven by input signals with r < 1, the analytical expressions 
obtained are too cumbersome to be of practical value. Thus, in the present work, we rely on 
the numerical treatment of the Langevin equation, Eq. (|2*|). following the procedure detailed 
in Ref. Q- 

In Fig.^ we depict the behavior of R ou t with the noise strength D, for input signals of the 
type given in Eq. (JTJ), with subthreshold amplitude A = 0.35, fundamental frequency Q = 
0.0024 and r = 0.1, 0.4, 0.7, 0.95, 0.98, 1. For all values of r, the characteristic nonmonotonic 
behavior of the SNR with D is obtained. As r increases, the maximum value of R out increases. 
Namely, the longer the potential remains asymmetric during each half-cycle, the larger the 
maximum height in the SNR is. Therefore, R out is quite sensitive to the duration of the 
pulses within each half-cycle. 
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FIG. 1: The output signal-to-noise ratio, R ou t vs. the noise strength, D, for r = 0.1 (circles), 
r = 0.4 (squares), r = 0.7 (triangles), r = 0.95 (filled circles), r = 0.98 (filled squares) and r = 1.0 
(filled triangles). The input signal has a subthreshold amplitude A = 0.35 and a fundamental 
frequency £1 = 2vr/T = 0.0024. 

In Fig. 121 the behavior of G with the noise strength D is depicted for the same values 
of r as in the previous figure. For all values of r considered, there exists a range of noise 
values such that G is larger than unity. The peak value of the SR gain has a nonmonotonic 
behavior with r. At the lowest value of r considered, the SR gain has a rather large peak 
value. Then, as r is increased, the peak of the SR gain decreases, in agreement with the 
observations in [3, Isj]. As the duration of the pulses gets larger so that r is closer to 1, the 
tendency of the SR gain maximum reverses and a considerable increase in the maximum is 
observed. 

The just mentioned features can be rationalized by noting that the SR gain depends 
on R out and on R in . As indicated above, for fixed values of the noise strength, D, and 
amplitude, A, the input Ri n always increases with r. Also, for given values of A and r, R in 
decreases monotonically with D. On the other hand, the results depicted in Fig. Q indicate 
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that there are two main effects on the location of the maximum of R ou t as r increases. First, 
as noted before, the maximum height increases as r increases. Second, the maxima appear 
at increasingly large values of D as the duration of the pulses increases. This second effect 
manifests itself clearly for pulses of sufficiently long duration, namely for r values larger 
than r m 0.9, while it is almost unnoticeable for smaller bvalues of r. For low values of r (let 
us say r = 0.1), even though the peak of R out is the smallest one appearing in Fig. [TJ the 
corresponding value of Ri n is so small (due to the smallness of r) that the SR gain reaches 
the large values depicted in Fig. |21 As r increases, the height of the R out maximum also 
increases, but appearing at an approximately constant value of the noise strength. Thus, the 
increase of Ri n with r counterbalances the increase of R out in such a way that the SR gain 
decreases. Finally, for long duration pulses, the shift to the right of the R out maximum and 
the large increase in its height are the cause of the increase in the maximum gain observed 
in Fig. El 

In conclusion, we observe that the behavior of the SR gain for pulses of relatively short 
duration is basically a consequence of the behavior of Ri n , in agreement with the observation 
of Gingl et al. [f], Q, Q|. On the other hand, for r close to 1, the behavior of the SR gain is 
dominated by the output SNR. 
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